Simulations of Helix Unwinding in Ferroelectric Liquid Crystals 
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In bulk ferroelectric liquid crystals, the molecular director twists in a helix. In narrow cells, this 
helix can be unwound by an applied electric field or by boundary efi'ects. To describe helix unwinding 
cis a function of both electric field and boundary effects, we develop a mesoscale simulation model 
based on a continuum free energy discretized on a two-dimensional lattice. In these simulations, we 
determine both the director profile across the cell and the net electrostatic polarization. By varying 
the cell size, we show how boundary effects shift the critical field for helix unwinding and lower the 
saturation polarization. Our results are consistent with experimental data. 
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I. INTRODUCTION 

One of the most extensively studied phases of liquid 
crystals, both for basic research and for applications, is 
the smectic-C* (SmC*) phase of chiral molecules. In this 
phase, the molecules lie in layers and are tilted with re- 
spect to the layer normal direction. The combination 
of molecular chirality, smectic layering, and tilt order 
leads to two effects: a ferroelectric polarization within 
the smectic layer plane and a helical modulation in the 
orientation of the molecular tilt from layer to layer y|. 
The ferroelectric polarization is useful for display devices, 
which use an applied electric field to switch the molecular 
orientation 01 ■ It is also useful for thermal sensors, which 
measure the temperature variation of the polarization, 
known as the pyroelectric effect 0, Q ■ Both of these ap- 
plications require a uniform orientation of the molecules. 
Hence, the helix must be suppressed, or unwound, by an 
applied electric field or by boundary effects. 

Helix unwinding has been modeled through continuum 
elastic theory. In a bulk SmC* phase, unwinding induced 
by an applied electric field can be described by the sine- 
Gordon equation, presented below d, IE 0- Under an 
electric field, the helix distorts and the helical pitch in- 
creases. If the field is increased above a critical thresh- 
old, the pitch diverges and the helix is suppressed. By 
contrast, helix unwinding induced by boundary effects 
in a narrow cell is more complex |E M ^E ^3 ■ ^^^^ 
case, a helix must have a series of disclination lines near 
the cell surfaces, which separate a twisted interior region 
from a uniform surface region. Continuum elastic the- 
ory shows that the unwinding transition is controlled by 
an energetic competition between the helical state with 
disclinations and the uniform untwisted state. If the cell 
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thickness is below a critical threshold comparable to the 
helical pitch, the uniform state is favored and the helix 
is suppressed. (A third possible mechanism for helix un- 
winding is shear flow jlJl, but this is not often used in 
SmC* liquid crystals, and we will not discuss it here.) 

In this paper, we investigate the unwinding of a SmC* 
helix by boundary effects, or by combined boundary and 
electric-field effects, through a series of Monte Carlo sim- 
ulations of a continuum free energy discretized on a 2D 
lattice. These simulations serve two purposes. First, 
they allow us to visualize the complex director configu- 
ration within a narrow cell as a function of cell thickness 
and electric field. We obtain snapshots of the molecular 
tilt profile through the helix unwinding process. Second, 
they allow us to relate the microscopic director config- 
uration to two macroscopic variables, an average chiral 
order parameter and the net electrostatic polarization of 
the cell. The latter variable can be compared with ex- 
perimental measurements of SmC* cells. 

In these simulations, we first consider boundary effects 
alone. We simulate the SmC* phase in narrow cells, and 
determine the structure of the helix as a function of cell 
thickness. We confirm that the helix unwinds at a thresh- 
old thickness approximately equal to the helical pitch, in 
agreement with continuum elastic theory. We then use 
the same Monte Carlo approach to simulate helix un- 
winding due to the combined effects of cell boundaries 
and electric field. In these simulations, we calculate both 
the molecular tilt profile across the cell and the net elec- 
trostatic polarization. We observe three distinct regimes 
of response, in which the helical pitch is first distorted 
and then expelled as the applied field is increased. By 
varying the cell thickness, we determine how boundary 
effects shift the critical field for helix unwinding and lower 
the saturation polarization. These simulation results are 
consistent with trends observed in experiments. 

The plan of this paper is as follows. In Sec. II, we 
review the continuum elastic theory of the SmC* phase. 
We show how an applied electric field unwinds the SmC* 
helix in a bulk system, and sketch an approximate en- 
ergetic argument for a finite cell. In Sec. HI, we use 
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FIG. 1: Idealized bookshelf geometry of a SmC* liquid crystal 
in a narrow cell. The helical axis is along the 2-axis, perpen- 
dicular to the smectic layers. An electric field is applied along 
the j/-axis. 

Monte Carlo simulations to investigate the effects of fi- 
nite cell thickness under zero applied electric field. In 
Sec. IV, we combine the effects of cell boundaries and 
applied electric field to determine the tilt profile and the 
polarization, and compare the polarization with experi- 
mental measurements. 



II. THEORY 

A SmC* liquid crystal in a narrow cell has the idealized 
bookshelf geometry shown in Fig. ^ The molecules lie 
in layers, and they are tilted with respect to the smectic 
layer normal. The 3D orientation is represented by the 
director n. This unit vector can be written in spherical 
coordinates as n = (sin 61 cos sin 6* sin 0, cos 61), where 
6 is the polar angle of the tilt and (j> is the azimuthal 
angle. The molecular tilt is conventionally described in 
terms of the projection of n into the layer plane, c = 
(sin 9 cos (j>, sin 6 sin cj)) . 

In this geometry, we expect the director to depend on 
the y and z coordinates. The z coordinate is along the 
smectic layer normal. Because of the molecular chirality, 
the molecular orientation rotates in a helix from layer to 
layer, which makes the director a periodic function of z. 
The y coordinate is the narrow dimension of the cell, with 
a thickness d of order microns, across which an electric 
field E is applied. The molecules may interact strongly 
with the front and back surfaces of the cell, at y — and 
d. As a result, the director may rotate as a function of 
y. By contrast, we do not expect the director to depend 
on the third coordinate x. In this bookshelf geometry, 
the system is uniform as a function of x. (In certain 
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FIG. 2: Unwinding of a SmC helix in a bulk system under an 
applied electric field E. Solid line: Perfect helix for E = 0. 
Dashed line: Distorted helix under moderate E. Dot-dashed 
line: Single soliton in the director, near the threshold for helix 
unwinding. Dotted line: Unwound uniform director. 

liquid cr ysta ls, the smectic layers buckle as a function 
of X [T^ Il4 ITslll^ . but we do not consider that effect 
here.) For that reason, we write the angles 9 and (j) as 
functions of y and z. 

We can now construct the simplest model free energy 
to describe variations in the director. This model must 
include four interactions. First, the smectic layer order 
interacts with the molecular orientation, and favors a 
particular tilt of the molecules with respect to the layers. 
This interaction can be expanded as a power series in the 
tilt magnitude |c| = sin0, which gives — ^rjcp + |u|c|^, 
for some series coefficients r and u. In terms of these 
coefficients, the favored tilt is |c| = {r/uY^^. Second, 
there is an electroclinic interaction of the molecules with 
the applied electric field E. Because of the molecular 
chirality, this interaction couples the field in the y di- 
rection with the tilt in the x direction, so it can be 
written as 6z • E x c = —bEyCx- In other words, the 
electrostatic polarization is P = &z x c, or Py — —bcx- 
Third, there is a chiral interaction that favors a variation 
of the director from layer to layer, which can be writ- 
ten as — Az • c X dc/dz. Fourth, there is the Frank free 
energy for elastic distortions of the director, which lim- 
its the variations from layer to layer. This contribution 
can be written as ^K{diCj){diCj), summed over i and 
j. Putting these pieces together, the total free energy 
density becomes 

If dc 1 

F = --r\c\^ + -u\c\'^+b7.-'Exc-\7.-cx—+-K{diCj){d,Cj). 

(1) 

This free energy is invariant under rotations in the xy 
plane, but it is not invariant under reflections, because 
reflection symmetry is broken by the molecular chirality. 

We can see immediately that this minimal model leads 
to a helix in the z direction. Consider a director of the 
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form c = (c cos 0, c sin 0) , where the magnitude c = sin 6 
is constant and the azimuthal angle depends on posi- 
tion. The free energy density simplifies to 

F = --rc^ + -uc'^ - hEccos(t> - Ac^^ + -Kc'^ |V<?!>|^ . 
2 4 az 2 

_ , (2) 

In the limit of zero electric field, i? = 0, we minimize this 
free energy over (f> to obtain 

^ - A 
'dz~K' 



(3) 



Hence, the azimuthal angle increases linearly as (jjiz) — 
4>Q + q^z, with go = A/if, as shown by the solid line 
in Fig. 121 This linear increase in (j) implies a perfect 
sinusoidal helix in c. 

This minimal model also predicts unwinding of the he- 
lix under an applied electric field. This behavior is anal- 
ogous to the theory of Meyer for a cholesteric phase in an 
electric or magnetic field [HSH- 
the free energy over </> gives 

dz"^ 



hE 
'k'c' 



For E ^ 0, minimizing 
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This is the standard sine-Gordon equation. The form 
of the solutions depends on the value of i?, as shown in 
Fig. 12 For low E, the helix is distorted, so that the di- 
rector is approximately aligned with the field in most of 
the system. As E increases, the helix becomes even more 
distorted, with sharper steps between domains where (p is 
approximately a multiple of 27r. Eventually the system 
crosses over into regime of uniform domains separated 
by sharp domain walls, or solitons. In that regime, a sin- 
gle soliton has the profile (j){z) = 2tt — 4tan~-'^[exp(— (z — 
•2waii)/C)], where ^ = (2Kc/bEy/^. As E continues to in- 
crease, the spacing between the domain walls increases, 
or equivalently their density decreases. At the critical 
threshold E = {n^ /8){X^c/bK), the last domain wall van- 
ishes and the system becomes uniform. 

The question is now: What other types of influences 
can also unwind a SmC* helix? Clearly one possibility 
is surface effects. Interactions along the front and back 
surfaces of a finite cell, at ?/ = and d, can anchor the 
director at those surfaces. If the elastic interactions de- 
scribed by the parameter K are sufficiently strong, and 
the cell thickness d is not too big, then this anchoring 
may extend throughout the interior of the cell, giving a 
uniform director. This is the basis of surface-stabilized 
ferroelectric liquid crystal cells 0]- 

The threshold thickness for unwinding a helix is not 
obvious. As shown in Fig. ^ the helical pitch is along 
the z direction, but the narrow dimension of the cell is 
along the y direction. Because these directions are per- 
pendicular, there is no simple geometric reason why a 
helix must unwind when the cell thickness is less than 
the pitch. Rather, there must be some energetic argu- 
ment that relates these two length scales. 

An energetic argument has been developed for narrow 
cells of the cholesteric phase and has been extended 



to narrow cells of the SmC* phase [lol[lll |. In its simplest 
form, the argument can be stated as follows . If a cell 
has a helix in the interior, but a uniform director along 
the front and back surfaces, then it must have a series of 
disclination lines running along the x direction near the 
surfaces. There must be one disclination line for each 
helical pitch. We can compare the energy of the helix 
(with disclinations) with the energy of the uniform state 
to find the threshold thickness for helix unwinding. The 
energy of the helix (with disclinations) is the negative 
energy gained from the helix plus the positive energy lost 
to the disclinations, 
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(5) 



where -Eunc is the disclination line energy per unit length. 
The helix unwinds if AE > 0, which implies 



d> 



Eli, 



K 



(pitch). 



(6) 



Since the line energy Eu^c should be of order K, the 
threshold thickness should be comparable to the pitch. 

In the following sections, we test this argument 
through a series of Monte Carlo simulations. In these 
simulations, we obtain snapshots of the director configu- 
ration for different cell thicknesses, both at zero field and 
under an applied electric field. These snapshots provide 
specific illustrations of the disclinations discussed above. 
For zero field, the simulations confirm that the helix un- 
winds at a critical thickness approximately equal to the 
helical pitch. For finite field, the simulations show helix 
unwinding induced by the combined effects of boundaries 
and electric field in a cell above the critical thickness. In 
both cases, we relate the microscopic snapshots of the 
director configuration to macroscopic variables. One of 
these variables, the net electrostatic polarization, can be 
compared with experimental measurements. 



III. FINITE CELLS UNDER ZERO ELECTRIC 
FIELD 

To model helix unwinding in a finite cell of the SmC* 
phase, we map the system onto a 2D square lattice. The 
lattice dimensions represent the yz plane shown in Fig. ^ 
y is the narrow dimension of the cell and z is the smectic 
layer normal. We assume the system is uniform in the x 
direction. On each lattice site (y, z) there is a 3D unit 
vector n(?/, z) representing the local molecular director. 
This vector can be parametrized in terms of the polar 
angle 0{y, z) and azimuthal angle ^(y, z), or equivalently 
in terms of the projection c(j/, z) into the smectic layer 
plane. 

For the lattice simulations, we discretize the free en- 
ergy of Eq. to obtain 



F = 



-ir|c(y,z)|2 



iu|c(2/,z)|* + 6z-E X c{y,z) 
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c(y,z) + c(y,z + l) , , ^, , 

-Az- '-x[c{y,z + l)-c{y,z)] (7) 



+ -i^(|c(y + 1, z) - c(y, + |c(y, z + 1) - c{y, z)p) 

with E = in this section. This free energy is similar to 
the free energy studied in Ref. 16] but with one impor- 
tant difference: that paper simulated the xy plane, but 
we now simulate the yz plane. 

A further consideration for the simulations is bound- 
ary conditions. Experimental cells may be symmetric, 
with the local director aligned along the same direction 
on both front and back confining walls, or they may be 
asymmetric, with the director aligned along one wall and 
an open boundary on the other side. In our simulations 
we use an aligning boundary condition with the director 
fixed on the wall at y = with a specified tilt angle. 
On the other wall at y = d we use the boundary con- 
dition 9c(j/, z)/dy = 0. This arrangement can represent 
an asymmetric cell, or one half of a symmetric cell, with 
the other half a mirror image of the first. Experimental 
cells are very large in the z direction so that the top and 
bottom boundaries should not affect the physics of the 
interior. In the simulations, we use the boundary condi- 
tion 9c(j/, z)ldz — for the top and bottom boundaries. 

We simulate the system with the parameters r = 
0.007596, u = 1, 6 = 1, £: = 0, if = 1.5, and A = 0.25 and 
0.125. The small value of r corresponds to a tilt angle of 
approximately 5°. We use the large z dimension of 160, 
and several values of the thickness d in the y direction. 
For each set of parameters and system size, we begin the 
simulations in a disordered state at a high Monte Carlo 
temperature, and then slowly reduce the temperature un- 
til the system settles into a single ground state and the 
fluctuations in c become negligible. This procedure can 
be regarded as a simulated-annealing minimization of the 
lattice free energy of Eq. ^ . 

To visualize the simulation results, we draw the yz 
plane in Figs. |31and2| The director configuration is rep- 
resented by short lines that show the projection of the 
3D director into the yz plane. Hence, vertical lines indi- 
cate Cy = 0, and tilted lines indicate Cy ^ 0. Because the 
lines representing the director are drawn close together, 
helical regions resemble twisted ribbons. 

Figure O shows the simulation results for A — 0.25 and 
K — 1.5. For these parameters, the favored wavevector 
is go = ^/K « 0.17, and hence the unperturbed pitch is 
w 38. For a thickness of 10, the system is uniform, 
with Cj, = everywhere except near the top and bot- 
tom surfaces. Those distortions are edge effects within a 
fractional pitch of the free surfaces, which do not affect 
the bulk behavior inside the cell. Hence, we see that the 
system of thickness 10 is unwound. By contrast, for a 
thickness of 12, the system shows a well-defined helix, 
with a periodic modulation of Cy throughout the cell, ex- 
cept very close to the aligning surface at y = 0. Thus, 
there is a clear helix winding/unwinding transition as a 
function of thickness. The transition occurs at a critical 
thickness between 10 and 12 for a cell with asymmetric 
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FIG. 3: Director configuration in the yz plane for simulations 
witli A = 0.25, leading to a pitch of approximately 38. The 
system has an aligning boundary along the left side (y = 0). 
Note that the y axis is exaggerated compared with the z axis. 

(a) For a thickness of 10 (y = through 9), the system is 
uniform, except for some edge effects at the top and bottom. 

(b) For a thickness of 12, the system has a clear helix. 
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FIG. 4: Director configuration in the yz plane for simula- 
tions with A = 0.125, leading to a pitch of approximately 75. 
(a) For a thickness of 23, the system is uniform, except for 
some edge eflects. (b) For a thickness of 25, there is a clear 
helix. 
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FIG. 5: The chiral order parameter x defined in Eq. ||HJ, 
as a function of the system thickness d. For A — 0.25 (pitch 
~ 38), the unwinding transition occurs at a thickness between 
10 and 12. For A — 0.125 (pitch « 75), the transition occurs 
at a thickness between 23 and 25. 



boundary conditions (or a half-thickness between 10 and 
12 if the simulation is regarded as half of a symmetric 
cell). This critical thickness is not equal to the pitch, 
but it is certainly the same order of magnitude, in agree- 
ment with the theoretical expectation of Eq. (jSJ. 

Figure^shows the corresponding results for A — 0.125. 
In this case, the favored wavevector is go ~ 0.083, so the 
unperturbed pitch is 27r/go ~ 75. For a thickness of 23, 
the system is uniform, again except for some edge effects 
near the top and bottom surfaces. By comparison, for 
a thickness of 25, there is a distinct helix throughout 
the interior of the cell. Thus, there is a helix wind- 
ing/unwinding transition with a critical thickness be- 
tween 23 and 25. This critical thickness is approximately 
twice the critical thickness of the previous case. Hence, 
we see that the critical thickness is approximately pro- 
portional to the pitch, again in agreement with the the- 
oretical expectation. 

For a quantitative measurement of helix winding and 
unwinding, we must define a chiral order parameter. One 
simple choice of a chiral order parameter is just the chiral 
term of the free energy © , without the factor of A itself. 
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X 



ciy,z)+c{y,z + l) 



X [c(y,z + l)-c(y,z)]. (8) 

Figure El shows this order parameter as a function of the 
system thickness d for both series of simulations, with 
A = 0.25 and 0.125. For A = 0.25, the plot shows a sharp 
transition between thickness 10 and 12, as x jumps from 
0.022 to 0.24. For A = 0.125, there is a distinct transition 
between thickness 23 and 25, as x jumps from 0.026 to 
0.13. This analysis confirms that the winding/unwinding 



transition occurs at a thickness that is proportional to, 
and the same order of magnitude as, the helical pitch. 



IV. FINITE CELLS UNDER AN ELECTRIC 
FIELD 



In the previous section, we showed that a SniC* helix 
can be unwound by surface effects in a finite cell. If 
the thickness is greater than the critical threshold, the 
helix is present at zero electric field. However, when an 
electric field is applied, the helix can be unwound by the 
combined effects of the surfaces and the electric field. In 
this section, we simulate that combination of surface and 
field effects. 

For these simulations, we use the same Monte Carlo ap- 
proach as in the previous section. We use the discretized 
free energy of Eq. (|HJ| with the parameters r = 0.0625, 
u = I, b = 1, K = 1.5, and A = 0.25. This value of 
r corresponds to a tilt angle of approximately 15°. We 
perform the simulations for four values of the cell thick- 
ness, d — 10, 20, 40, and 60, and scan through the elec- 
tric field E at each thickness. For these parameters, the 
unperturbed pitch is 2t:K/X « 38 and hence, based on 
the results of the previous section, the zero-field wind- 
ing/unwinding transition occurs at a thickness between 
10 and 12. Hence, the simulations at thickness 10 should 
be unwound at all values of the electric field, while the 
simulations at larger thickness should go through the 
winding/unwinding transition as a function of electric 
field. 

Figure IHl shows the simulation results for the system of 
thickness 20. For E = 0, the system has a helix every- 
where in the cell, except a narrow region near the aligning 
boundary at y = 0. By comparison, for E = 0.004, the 
helix is suppressed in much of the cell. It persists only in 
regions of the cell near the free surfaces at the top and 
bottom. This behavior near the free surfaces is consis- 
tent with the previous section, which showed that free 
surfaces tend to favor the helical modulation. When the 
field increases to E — 0.006, the helix is suppressed in 
more of the cell, and it persists only in smaller regions 
near the top and bottom surfaces. Once the field reaches 
E = 0.008, the helix is suppressed throughout the inte- 
rior of the cell. The director in the cell is now uniform, 
except for very narrow edge effects at the top and bot- 
tom. Hence, the electric field has driven the finite cell 
through the helix winding/unwinding transition. 

As in the previous section, we need an order parameter 
to describe the wiirding/unwinding transition quantita- 
tively. In this case, the electrostatic polarization provides 
an experimentally relevant order parameter, which shows 
how the net polar order of the cell couples to the elec- 
tric field. As argued in Sec. II, the polarization is the 
quantity conjugate to the electric field in the free energy, 
and hence P = &z x c, or P,, 



-bcx- We average this 
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FIG. 6: Director configuration in the yz plane as a function of increasing electric field, for a fixed thickness of 20. (a) E — 0. 
(b) E = 0.004. (c) E = 0.006. (d) E = 0.008. 



quantity over the cell to obtain 

-''sites - . 

Figure [7\ shows the simulation results for the polariza- 
tion as a function of electric field for each cell thickness. 
Figure Et^) presents the results on a linear scale, and 
Fig.[7fb) presents the same results on a logarithmic scale. 
Note that the polarization is not zero at zero field because 
of the symmetry-breaking surface alignment at y = 0. 

From the plots in Fig. [7| we can see that the system 
has three distinct regimes in its response to an electric 
field. For low field E ^ 0.01, there is a regime of he- 
lix unwinding. In this regime, an increasing electric field 
gradually aligns the directors, suppresses the helix, and 
prevents the local polarization from averaging to zero. As 
a result, the net polarization increases rapidly as a func- 
tion of field. (This low-field regime does not occur for the 
lowest thickness d = 10, because the helix is suppressed 
by surface effects even without a field.) For intermedi- 
ate field 0.01 ^ E < 1, there is a regime of suppressed 
twist. In this regime, the helix is already unwound, so the 
only effect of electric field is to increase the local tilt and 
polarization. Hence, the net polarization increases more 



slowly, roug hly as E^/^. Finally, for high field E > 1, 
there is a saturated regime. Here, the helix is already 
unwound, the local tilt is at its maximum value Cx = 1, 
and the polarization is at its maximum value of Py — b. 
Although the simulations only go to £' = 1.2, we see that 
the polarization cannot increase at higher field because 
it is already saturated. 

In Fig. we can also compare the relative polariza- 
tion of thinner and thicker cells. For low field, thinner 
cells have a higher polarization than thicker cells, because 
the helix reduces the polarization in thicker cells but the 
aligning boundary at ?/ = unwinds the helix in thin- 
ner cells. By contrast, for high field, thicker cells have 
a higher polarization than thinner cells, because there 
is no helix at any thickness, and surface effects at y = d 
suppress the polarization in thinner cells. A similar high- 
field limit has been discussed by Shenoy et al. ^T§\, who 
see the same effect experimentally in cells with different 
boundary conditions. For intermediate field, the polar- 
ization of thinner and thicker cells must cross. 

To compare our results with experimental measure- 
ments of the polarization as a function of electric field, 
we must take into account one subtlety of the experi- 
ments. As shown by Ruth et al. 0,1131, the polarization 
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FIG. 7: Simulation results for the polarization as a function 
of electric field in cells of four thicknesses, (a) Linear scale, 
(b) Logarithmic scale. On the logarithmic scale, note the 
three regimes of helix unwinding, suppressed twist, and satu- 
ration. 

observed experimentally (by the triangle-wave technique 
or other techniques) is not the total polarization conju- 
gate to the electric field. Rather, it is a specific nonlinear 
component of the polarization, which can be written as 

PME) - P{E) - E . (10) 

The difference between the total polarization P{E) and 
the observable polarization Pobs(-E') is small when P{E) is 
saturated, but it is significant whenever P{E) varies with 
E. Hence, we must extract Pohs{E) from the simulations 
and compare that quantity with experiments. 

To extract Pohs{E) from the simulations, we need to 
calculate the derivative dP/dE. For that reason, we fit 
the simulation results for the polarization to the function 
P{E) = (1 -I- aE)/{(3 + jE). This is just an empirical 
fitting function, with no theoretical basis, but it gives 
a fairly good fit to the data in Fig. d^a). We can then 
differentiate this function and calculate the observable 
polarization Pohs{E) that corresponds to the simulation 
results. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Electric field (arb. units) 

400 T 




A , , , 1 

12 3 4 

Electric field (V/|im) 



FIG. 8; (a) Simulation results for the observable polarization 
fobs ~ P — E{dP/dE) in cells of four thicknesses, (b) Exper- 
imental data for Pobs in three cells I2H . 

Figure ISJa) shows the results for Pobs {E) extracted 
from the simulations at each system size. This function 
is small for low E, then increases towards its saturation 
level Pobs = b at high E. As discussed above, thinner 
cells have a higher value of Pobs at low field, and thicker 
cells have a higher value of Pobs at high field. By compar- 
ison, Fig.|SJb) shows sample experimental measurements 
of Pobs as a function of applied electric field in three cells 
with asymmetric boundary conditions |2l| . The material 
is 10PPBN4 (described in Ref. HI), and the tempera- 
ture is 7 °C below the SmA-SmC transition. Note that 
the experimental data show the same general features 
as the theoretical plots. In particular, we see the same 
crossover between higher Pobs in thinner cells at low field 
and higher Pobs in thicker cells at high field. Thus, the 
simulation results are consistent with the trends seen in 
experiments. 

In conclusion, we have developed an approach for sim- 
ulating the helix winding/unwinding transition in SmC* 
liquid crystals. This approach is based on a minimal 
model for the free energy, which includes a chiral term 
that favors a helical modulation of the director from layer 
to layer. This bulk free energy competes with surface ef- 
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fects and clcctric-ficld effects, which both favor a uniform 
aUgnmcnt of the director. In zero field, the competition 
between the bulk helix and the surface alignment leads 
to helix unwinding at a critical thickness approximately 
equal to the helical pitch. When an electric field is ap- 
plied, the field-induced alignment adds to the surface ef- 
fects, and induces helix unwinding even for thicker cells. 
The electrostatic polarization is an appropriate order pa- 
rameter to quantify this field-induced unwinding, and our 
simulation results for the polarization are consistent with 
experimental measurements. 



Acknowledgments 



We would like to thank R. B. Meyer for helpful dis- 
cussions and K. A. Crandall for sharing his unpublished 
data with us. This work was supported by the U. S. Navy 
through Contract No. N00173-99-1-G015, and by the 
National Science Foundation through Grant No. DMR- 
9702234 and DMR-01 16090. 



[1] P. G. de Gennes and J. Prost, The Physics of Liquid 
Crystals, 2nd edition (Oxford, Oxford University Press, 
1993). 

[2] J. W. Goodby, R. Blinc, N. A. Clark, S. T. Lagerwall, M. 
A. Osipov, S. A. Pikin, T. Sakurai, K. Yoshino, and B. 
Zeks, Ferroelectric Liquid Crystals: Principles, Proper- 
ties and Applications (Gordon and Breach, Amsterdam, 
1991). 

[3] J. W. O'SuIlivan, Yu. P. Panarin, and J. K. Vij, J. Appl. 

Phys. 77, 1201 (1995). 
[4] K. Crandall, D. Shenoy, S. Gray, J. Naciri, and R. 

Shashidhar, Proc. IRIA-MSS on Infrared Materials, p. 

173 (1999). 

[5] R. B. Meyer, Appl. Phys. Lett. 12, 281 (1968) 
[6] R. B. Meyer, Appl. Phys. Lett. 14, 208 (1969). 
[7] R. D. Kamion and J. V. Selinger, J. Phys.: Condens. 

Matter 13, Rl (2001). 
[8] P. E. Cladis and M. Kleman, Mol. Cryst. Liq. Cryst. 16, 

1 (1972). 

[9] M. Luban, D. Mukamel, and S. Shtrikman, Phys. Rev. A 
10, 360 (1974). 

[10] M. Glogarova, L. Lcjcck, J. Pavel, V. Janovcc, and J. 

Fousck, Mol. Cryst. Liq. Cryst. 91, 309 (1983). 
[11] T. Povse, I. Musevic, B. Zeks, and R. Blinc, Liq. Cryst. 



14, 1587 (1993). 
[12] A. D. Rcy Phys. Rev. E 53, 4198 (1996). 
[13] G. P. Crawford, R. E. Goer, J. Naciri, R. Shashidhar, 

and B. R. Ratna, Appl. Phys. Lett. 65, 2937 (1994). 
[14] R. E. Geer, S. J. Singer, J. V. Selinger, B. R. Ratna, and 

R. Shashidhar, Phys. Rev. E 57, 3059 (1998). 
[15] F. J. BartoU, J. R. Lindle, S. R. Flom, R. Shashidhar, 

G. Rubin, J. V. Selinger, and B. R. Ratna, Phys. Rev. E 

58, 5990 (1998). 
[16] J. V. Schngor, J. Xu, R. L. B. Selinger, B. R. Ratna, and 

R. Shashidhar, Phys. Rev. E 62, 666 (2000). 
[17] R. B. Meyer, private communication. 
[18] D. Shenoy, A. Lavarello, J. Naciri, and R. Shashidhar, 

Liq. Cryst. 28, 841 (2001). 
[19] J. Ruth, J. V. Selinger, and R. Shashidhar, Appl. Phys. 

Lett. 65, 1590 (1994). 
[20] D. Hermann, L. Komitov, S. T. Lagerwall, G. Anders- 
son, R. Shashidhar, J. V. Selinger, and F. Giesselmann, 

Ferroelectrics 181, 371 (1996). 
[21] K. A. Crandall, private communication. 
[22] J. Naciri, B. R. Ratna, S. Baral-Tosh, P. Keller, and R. 

Shashidhar, Macromolecules 28, 5274 (1995). 



